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Abstract: In this paper we consider a special kind of three-point functions of HHL type 
at weak coupling in AA = 4 SYM theory and analyze its volume dependence. At strong 
coupling this kind of three-point functions were studied recently by Bajnok, Janik and 
Wereszczynski [1]. The authors considered some cases of HHL correlator in the su(2) sector 
and, relying on their explicit results, formulated a conjecture about the form of the volume 
dependence of the symmetric HHL structure constant to be valid at any coupling up to 
wrapping corrections. In order to test this hypothesis we considered the HHL correlator in 
su(2) sector at weak coupling and directly showed that, up to one loop, the finite volume 
dependence has exactly the form proposed in [1]. Another side of the conjecture suggests 
that computation of the symmetric structure constant is equivalent to computing the cor¬ 
responding set of infinite volume form factors, which can be extracted as the coefficients of 
finite volume expansion. In this sense, extracting appropriate coefficients from our result 
gives a prediction for the corresponding infinite volume form factors. 
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1 Introduction 


It’s strongly believed that N = A SYM theory is integrable. There are a huge amount of 
supportive evidences in the literature starting from [2], where the relation between spin 
chain Hamiltonian and the dilatation operator of the theory was established, for detailed 
review, see [3]. Among one of the most exciting achievements revealing the power of the 
integrability is the so-called Quantum Spectral Curve (QSC) approach [4, 5] which provides 
efficient technique for computing anomalous dimension of a gauge invariant operators of the 
theory. 

Due to conformal symmetry, the computation of general correlation functions of AA = 4 
SYM can be reduced to computation of two- and three-point functions by means of the 
operator product expansion (OPE). In this sense these two quantities are the fundamental 
blocks of the theory. Two-point functions of gauge invariant operators, again due to the 
conformal symmetry, are completely dehned by the value of the corresponding anomalous 
dimension, and thus can be computed in terms of QSC technique. However for the moment 
there is no such an efficient analogue for computing three-point functions. Nevertheless 
signihcant progress have been done in this direction starting from [6, 7] and later developed 
systematically in [8], where the three-point functions of the su(2) sector at tree level were 
considered and, using the mapping between spin chains and gauge invariant operators, the 
computation of the structure constants was reduced to computation of the scalar products 
between two off-shell Bethe-states. This method was improved in [9] for some special 
conhguration, which allows to express the structure constant in terms of on-shell/off-shell 
scalar products and can be computed by Slavnov determinants. Later this result was 
extended to one loop [10, 11] and other rank-one sectors also were investigated [12, 13]. 
The great advancement was made very recently in [14], where the authors proposed all-loop 
procedure for the structure constant in AA = 4 SYM. Another interesting direction inspired 
by the light-cone string held theory was initiated in [15], where the authors interpreted 
the OPE coefficients as “generalized Neumann coefficients” and proposed a set of bootstrap 
axioms for these generalized Neumann coefficients. In parallel, the spin vertex approach 
was developed in [16-18] which can be seen as a weak coupling counterpart of the string 
vertex. 

Apart from direct approach, there is another method for the computation of three-point 
functions by means of their relations with the form factors. This approach was initiated 
in [19, 20], where the set of the axioms for the world-sheet form factor of the light-cone 
gauge fixed AdS^ x string theory was proposed. Recently in [1] the symmetric HHL 
(heavy-heavy-light) correlator at strong coupling of the su(2) sector was considered. The 
proposals for computing the HHL three-point functions were first formulated in [21, 22]. 
The authors of [1] showed that the prescription in [21] was inadequate and proposed an 
improved prescription. Using this new prescription the authors computed several examples 
of the three-point correlators and showed that their volume dependence exactly coincide 
with the hnite volume structure of the appropriate form factor. Guided by this result they 
proposed a conjecture that this hnite volume dependence should hold at any coupling of 
the theory. 
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The conjecture consists two parts. First, the finite volume dependence (neglecting 
wrapping) of a symmetric HHL structure constant is completely encoded into diagonal 
minors of the Gaudin determinant of a heavy state and have the following form dictated 
by results of form factor theory ([23]): 


Chhl 


/^ + EPiv(W)/^(i,-,t-,^) + EPiv({i,j})/®(i,..., 
_*__ 

Piv({l,...,iV}) 


i,...,iv) + ... 


( 1 . 1 ) 


and /9Ar({i, J, •••})’s are diagonal minors of the determinant: 

/9Ar({l,---,A^}) = dejt f ^ (^Lpfc(u) + ^ - log ^(riz, u^)) j. (1.2) 

Second part of the conjecture, also coming from the form factor approach, is that the 
coefficients in the expansion (1.1) are diagonal infinite volume form factors of a light 
operator O. In this paper we test the first part of the conjecture at weak coupling limit and 
show that the form (1.1) of the finite volume dependence holds at tree level as well as at one 
loop. At the same time the second part of the conjecture allowed us to predict expressions 
for diagonal infinite volume form factors by extracting corresponding coefficients from the 
expansion (1.1). 

The structure of the paper is the following. In section 2 we set up the problem. In 
section 3, we will describe in detail the map between the field theory operator Oa and the 
spin chain operator Oq,. In section 4, we review the integrability tools that we need, which 
includes algebraic Bethe ansatz (ABA), the scalar products between Bethe states and the 
solution of quantum inverse scattering problem (QISP). In section 5, we review the form 
factor formalism and the finite volume corrections. In section 6 and 7, we give a proof of 
the finite volume structure of the three-point function conjectured in [1] at tree level. In 
section 8 we compute some examples of the symmetric HHL correlator and list expressions 
for the corresponding infinite volume form factors. In section 9 we give general arguments 
proving that the same finite volume structure holds as well at one loop. Finally we conclude 
and discuss future directions to pursue in section 10. 


2 Set-up 

The central object of our paper is the symmetric HHL three-point correlation function 
of gauge invariant operators, by which we mean that two heavy ^ (the length L ^ 1) 
operators are conjugated to each other and the third operator has a few number of sites. 
The technique, which we are going to use, is due to the relation between the gauge invariant 
operators and the spin chains. Following [24], the symmetric structure constant can be 
mapped to diagonal matrix element of a spin chain operator. We take the heavy operators 
to be the eigenstates of the su(2) sector dilatation operator and conjugated to each other. 
They are constructed from the following scalar fields 

^Although we call it "heavy" the results obtained in this paper are valid for small values of L as well 
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Oi : {Z,X} 




(2.1) 


At tree-level, we consider operators with definite one-loop anomalous dimension [8]. The 
one-loop dilatation operator in the su(2) sector is the Heisenberg XXX 1/2 spin chain Hamil¬ 
tonian whose eigenstates can be constructed by Bethe Ansatz techniques. Since we are con¬ 
sidering diagonal matrix elements, the wave functions of the two operators are conjugate 
to each other. The third operator will be denoted hereafter by Oa- In this paper, we con¬ 
sider Oa in the compact sector. Only the zero i?-charge terms give non-zero contribution. 
Therefore, we are interested in the operators of the following form 

TrZZ, TrXX, TrXZZX, ••• (2.2) 


This kind of operators are in the 5o(4) sector of A/" = 4 SYM theory. In addition, we require 
Oa have definite anomalous dimension. The three-point function is fixed by conformal 
symmetry up to the structure constant Ca 


(Ol(xi)Oa(x2)Ol(x3)} 


La A/i\/A/’q,C'q, 

Nc kl2|^12|a;^3|Ai3|a;23|A23 


(2.3) 


where 


u u 

^ij = 


— X- 


Aij = -(Ai + Aj - Ak 


(2.4) 


and L is the length of Oi while La is the length of the operator Oa- The two-point functions 
are normalized as 

(O,(x0O,(x2)) = (0„(I1)0„(X2)> = (2,5) 

\Xl2\ ^ \Xl2\ 

The structure constant Ca can be expressed in terms of correlation functions of the Heisen¬ 
berg spin chain 


_ (u|d„|u) 
(u|u) 


( 2 . 6 ) 


where |u) = - denotes an on-shell Bethe state, corresponding to the heavy 

operator, and Oa = Oa{crf ^af) is an operator made of local spin operators. In this way, 
the computation of three-point function in planar N = A SYM theory is recast into the 
calculation of correlation functions in the Heisenberg spin chain. 


3 Prom field theory correlation functions to spin chain matrix elements 


In this section, we summarize how to write the field theoretic operators Oa in terms of spin 
chain operators Oa- Let us introduce the following notation 


z = 4P, x = 4>\ 


z = 
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(3.1) 






The light operator is the linear combination of the single trace operators Tr cj)^^ ■ ■ ■. 
By planarity, only operators of the following form will contribute to the three-point function 

Tr(/)*/•• •(/>*'4,jfc=0,1. (3.2) 

where the indices 1, • • • ,21 denotes the position on the third spin chain and 21 = La is the 
length of the third operator. The zero i?-charge condition is given by 

l 

^{in-jn)=0. (3.3) 

n=l 

It is not hard to see that the operator (3.2) can be mapped to the following spin operator 

|24] 

where the indices n, • • • ,n — I + 1 denote the positions on the long spin chain, and the 
operators are the basis 2x2 matrices (E“^)jj = Sai^bj in the local quantum space 
LLn = C^. The operators E“^ are related to the local spin operators as follows 

= + e12^ct+, E22^1(I-a^). (3.5) 

Here a^, are the usual Pauli matrices acting on the space spanned by If) and IJ.). By 
the mapping (3.4), we can translate the held theory operators into the spin operators. As 
an example, we consider the Konishi operator 

Ok = Tj:XX+ Ti:YY+ TiZZ. (3.6) 

Since the heavy operators are in the su(2) sector, the contraction with Y and Y are zero 
and can be neglected. The Konishi operator can be mapped to the following spin operator 

Ok ^ Ok = E“ + Ef = L (3.7) 


Therefore we see that at tree level the structure constant, with the light operator being the 
Konishi operator, is trivial 


Ck 


(u|Oj^|u) 

(u|u) 


(3.8) 


In order to obtain non-trivial structure constant, the light operator must have at least 
La = 4. This corresponds to the insertion of two spin operators between the Bethe states. 
An example for length-4 operator is given in Fig. 3. 

The main focus of the current paper is the length-4 held theory operators, which cor¬ 
responds to operators in the spin chain which acts on two neighboring sites. These are the 
simplest non-trivial cases which can be studied thoroughly. We use and Ig to denote the 
lengths of held theory operators and spin chain operators, it is obvious that La = 2ls. 
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zzzxzxxzzzzxz 



Figure 1. An example of the mapping between field operators to spin chain operators. The 
operator in the field theory O = TrZAZX is mapped to the spin chain operator O = 

In order to have non-zero diagonal matrix element in (2.6), the spin chain operator 
shonld not alter the total spin of the state it acts on. At = 1, there are two independent 
operator, satisfying this condition, the identity and 

oi{n)=El^. (3.9) 

For Ig = 2, there are six independent operators: the nnity, oi(n), oi(n -|- 1) ^ and the other 
three operators are the following 

olin) = Ei^El\,. oUn) = EfE}^.,,. (3.10) 

For later convenience we also introdnce the following operator 

ol = E^X+i = (I - E“)(I[ - Eii+i) = I - oi(n) - oi(n + 1) + o^(n). (3.11) 

In what follows, we will stndy in detail the diagonal matrix elements of the spin operators 
Ol and O 2 (i = 1,2,3,4). 

4 The Algebraic Bethe Ansatz approach 

In this section, we review the algebraic Bethe Ansatz (ABA) approach to the correlation 
fnnctions in the XXX 1/2 Heisenberg spin chain, see [25-27] and references therein. This 
approach is based on two main elements: the Slavnov determinant formnla and the solution 

^Due to translation invariance of the spin chain at the level of correlation functions these two operators 
obviously coincide. 
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of the Quantum Inverse Scattering Problem (QISP). The Slavnov formula states that the 
scalar product of any Bethe state with an on-shell one can be written in terms of a deter¬ 
minant. The solution of the QISP enable us to express the local spin operators in terms of 
the elements of the monodromy matrix, which are non-local operators acting on the spin 
chain. 


4.1 Algebraic Bethe Ansatz 


In this subsection we give a brief review of the algebraic Bethe Ansatz method which also 
serves to fix our conventions. For more detailed and pedagogical discussions we refer the 
readers to [28]. The central object in ABA is the quantum i2-matrix. For the Heisenberg 
XXXi /2 spin chain, the ii-matrix is given by 


— U lafe T i Pafe 


^u-|-i00 0 ^ 

0 u i 0 
0 i u Q 
y 0 Q ^ u + ij 


ab 


(4.1) 


where lab Pab are identity and permutation operators acting on the space (8) C^, 
respectively. It has a complex parameter which is usually referred to as the spectral pa¬ 
rameter. Consider a spin chain of length L. At each site of the spin chain we define the 
Lax matrix 


Laniu dn) — Ran(u 9n ^/2)) — 1; ' ' ' j T (d'2) 

where n and a denote the quantum space and auxiliary space, respectively. The parameter 
9n is called the inhomogeneity associated to the site n. In the case of the homogeneous 
Heisenberg spin chain all the inhomogeneities are set to be zero, 9n = 0. The Lax operator 
obeys an important property that it satisfies the following braiding relation, called the RLL 
relation 


Rab(,^ '^)Tan(u)L5^('u) — Llyii(^V^Lan(^‘^^Rabi,'^ 


(4.3) 


Taking the ordered product of Lax operators at all sites, we obtain the monodromy matrix 

L 

Ta{u) = llLan{u-9n). (4.4) 

71=1 

It can be represented in the auxiliary space as a 2 x 2 matrix 


T{u) 


( A{u) B{u) \ 
\C{u) D{u) j • 


(4.5) 


By construction, the monodromy matrix is a non-local object and its elements are operators 
which act on the whole spin chain. From the RLL relation (4.3) and the definition of the 
monodromy matrix (4.4), one can show that the monodromy matrix satisfies a similar 
relation, called the RTT relation 


Rah{u - v)Ta{u)Tb{v) = Tb{v)Ta{u)Rab{u - v). 
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(4.6) 




Once written in terms of the components, the RTT relation leads to a quadratic algebra, 
between the ABA operators A, B, C, D. The most relevant relations for us are 

A{v)B{u) = f{u — v)B{u)A{v) + g{v — u)B{v)A{u), (4.7) 

D{v)B{u) = f{v — u)B{u)D{v) + g{u — v)B{v)D{u), 

[C{v), B{u)] = g{u — v)[A{v)D{u) — A{u)D{v)] 

where the functions f[u) and g{u) are 

.. XL , i ,, 

f(u) = -, g{u) = (4.8) 

u u 

All the conserved charges of the system are encoded into the transfer matrix, which is 
defined to be the trace of monodromy matrix in the auxiliary space 

Tiu) = Tra T{u) = A{u) + D{u). (4.9) 

The eigenstates of the transfer matrix simultaneously diagonalize all the conserved charges, 
including the Hamiltonian. In order to diagonalize it, one starts with a ferromagnetic 
reference state, called the psuedovacuum |H) = | f^), such that 

A{u)\n) = a(n)|H), D{u)\n) = d{u)\n), C{u)\n) = 0. (4.10) 

For later convenience, we define the Baxter polynomials 

L 

Qe{u) = W{u-On), Qg = Qe{u±i/2). (4.11) 

n=l 

Then we have 

a{u) = Qg{u), d{u) = Qg{u). (4-12) 

The eigenstates of the transfer matrix are generated by the action of a sequence of B- 
operators on the pseudovacuum 


u) = B{ui) • • •H(ri 7 v)|H) 


(4.13) 


provided that the rapidities u = {ui, ■ ■ ■ ,uisf} satisfy the Bethe Ansatz equation (BAE) 


ajuk) _ 1^ Uk -Uj +i 

d{uk) 

Uk 


(4.14) 


These elementary excitations are commonly referred to as magnons. If the rapidities satisfy 
(4.14), the corresponding Bethe state |u) is called on-shell, otherwise it is called off-shell. 
The on-shell Bethe states diagonalize the transfer matrix 


T(u)|u) = tu('w)|u) 
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(4.15) 





with the eigenvalue 


, f \ ( \Qu{u i) ,/ NQu('ti + 0 

tu{u) = a{u) + d{u)————. 

Qu{u) Qu[u) 

Here, similarly we have defined the Baxter polynomial for the rapidities 

N 

Qu{u )= 

k=l 


(4.16) 


(4.17) 


In this paper, we will compute the matrix elements of spin chain operators between two 
Bethe states. To this end, it is important to have a manageable expression for the scalar 
products between Bethe states. 

4.2 Slavnov determinant and Gaudin norm 

The scalar product of Bethe states is defined as 

N N 


(v|u) = (H| C{vj) B{uk)\n). 

j=i k=i 


(4.18) 


The Slavnov theorem [29] states that, if one of the Bethe states, say |u) is on-shell, the 
scalar product (4.18) is given by 


N 


(v|u) = a{vj)d{uj)yu,v, 
j=i 


where 


^u,v = 


detjk Q.{uj,Vk) 

uj-Vk+i 


The matrix element Q{uj,Vk) is given by 

Q{Uj,Vk) = t{Uj - Vk) - - Uj) 

. {uj -Vk-i) - (uj -Vk + i)e 


(4.19) 


(4.20) 


(4.21) 


2ipu ('«*;) 


= z 


{Uj - Vk)[iUj - VkY + 1] 

where Pu('i^) is the pseudomomentum 

2ip^(u) ^ d{u) Q^{u + i) 


a{u) Quiu - i)' 


and t{u) is 


/ N 1 1 

= T - 


(4.22) 


(4.23) 


u u + i 

In terms of the pseudomomentum, the Bethe Ansatz equation takes a particularly simple 
form 


^2ipu{uk) _ 
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(4.24) 









In the Slavnov formula, the rapidities {v} can be any set of complex numbers. Let us 
now consider the case when some of the Vk’s coincide with the Bethe roots, namely Vk = Uk 
for some k. The matrix element reads as 


Q{uj,Uk) 


2i 

{Uj - + 1 


/ TV 

I ^ ^ '7~~ 


2 

ui)'^ + 1 


L 

E 




(4.25) 


When all {v} coincide with {u}, the scalar product corresponds to the norm of the 
Bethe state. This norm is given by the Gaudin formula which takes the form of a Jacobian 
determinant. Let us define the following quantity 


N 

4'fc(u) = p{uk)L - i'^log S{uk -Uj) (4.26) 

where the momentum and the two-body scattering matrix of the Heisenberg spin chain are 


p{u) 



a{u) 
d{u) ’ 


S{u — v) 


u — V — i 
u — V + i 


In terms of the function <Lfc(u), the Bethe Ansatz equation (4.14) reads as 


(4.27) 


d>A:(u) = 27 r 4 , A; = l,...,iV 


(4.28) 


where G Z are the Bethe quantum numbers. The Jacobian for the change of variables 
between Uk and Ik is given by 


PAr({l,...,iV}) 


g^fc(u) 

duj 


det i Uk)- 
jk 


(4.29) 


The norm of an on-shell Wmagnon Bethe states is proportional to /3Ar({l, • • • , A^}), explic¬ 
itly 


(u|u) = j n a{uj)d{uj) JJ ^ ^ y j /9Ar({l,...,iV}) (4.30) 

\j=l j<k ^ ^ ^ J 

We will see below that the Jacobian /9jv({l; • • • j dl}) and its sub-diagonal counterparts take 
into account all the finite volume corrections of the diagonal matrix elements. 


4.3 The solution of quantum inverse scattering problem 

In the ABA approach, the idea to compute the diagonal matrix element (u|On|u) is to 
act the spin chain operator On on the bra state, so that the ket Bethe state |u) is left 
on-shell. The operator On is typically a multilocal operator in the sense that it acts on a 
finite interval of the spin chain. In order to apply the Slavnov formula, the bra state also 
needs to be a Bethe state, although we do not require it to be on-shell. The solution of the 
quantum inverse scattering problem (QISP) relates the local spin operators to the matrix 
elements of the monodromy matrix. These latter are non-local operators. In addition, from 
the Yangian algebra (4.7), it is clear that the action of A, B, C, D operators on a Bethe 


10 














state always gives a sum over Bethe states. Therefore, using the solution of QISP we can 
translate any local operator On into a sequence of ABA operators, and the state (u|0„ can 
be written as a sum of dual Bethe states. This enable us to apply the Slavnov formula 
discussed in the last subsection. We present the main statement of the solution of QISP in 
this subsection and refer to [26] for the proof and details. 

Let (a, 6 = 1,2) be the operators defined in (3.5) which act on the local quantum 
space Tin = C^. They can be represented by the elements of monodromy matrix as 


where 


' n—l 


=< n +*/ 2 ) \ n +'/ 2 ) 


.fc=l 


.k=l 


(4.31) 


T^\u) = A{u), T^^{u) = B{u), T^\u) = C{u), T^^{u) = D{u). (4.32) 


Once inserted inside a correlator, the transfer matrices in (4.31) act on a Bethe state and 
can be replaced by their eigenvalues 


(u|E“V> = 


(u|r“'’(0„ + i/2)|u) 


t\i{0n + */ 2 ) 

where t^x{u) is given by (4.16). Ai u = On + we have 

t\i{0n + i/‘^) — Qe{(^n + ^ 


QiiOr, 


(4.33) 


(4.34) 


The generalization of (4.33) to a string oi I + 1 operators is straightforward 


(u|E“o'o...e3|u)= m 


Qti^n+k) 


^_Q Q6{0n-\-k '^)Q\: 


(^n+Zc) 


Hn^n5(^n+fc+i/2)|u) (4.35) 


k=0 


We can now compute the r.h.s. of (4.35) by using the Yangian algebra and the Slavnov 
formula. 


5 Finite volume diagonal form factors 

The interest in the volume dependence of diagonal matrix elements of local operators is 
due to their appearance in different physical quantities. They are central objects in the 
calculation of finite temperature one-point functions [30, 31], as well as they play an essential 
role in the form factor perturbation theory [32]. The short distance behavior of correlation 
functions also involves the vacuum expectation values as basic ingredients [33]. Recently, 
it was conjectured [1] that the Heavy-Heavy-Light symmetric structure constants of the 
AdS/CFT correspondence (2.6) are also related to these quantities. This latter case is of 
the central interest of the current paper. To better understand this conjecture, in this 
section we summarize the theory of form factors in integrable models, starting from the 
infinite volume description and then take into account the finite volume corrections, up to 
wrapping. 
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5.1 Form factors in infinite volume 

We consider a 1 + 1 dimensional integrable quantum field theory defined by its S-matrix. 
For simplicity we restrict ourselves to diagonally scattering theories with a self-conjugated 
particle, the generalization to any diagonally scattering theory is straightforward. For a 
detailed review of the general case see [34], 

The inhnite volume states can be characterized by the set of momenta of particles. 
In 1 -|- 1 dimension it is convenient to label the particles by their rapidities m, the energy 
e{u) and momentum p{u) being a single valued functions. In the remote past, an in state 
consists of particles of ordered rapidities: the fastest one is the leftmost while the slowest 
is the rightmost. Contrarily, the particles in an out state are reversely ordered. 


ui, - ■ ■ , un) 


tti, • • • 

Ui,--- ,nAr)°"* 


Ui > ■ ■ ■ > UN 
ui < ■ ■ ■ < Un 


(5.1) 


The inhnite volume states that differ only in the order of rapidities are related by the two 
particle S-matrix^ 


I'i^l) ■ ■ ■ ) ) ^i+l) ■ ■ ■ ) ) — <5('Uj, Uj-i-i) I Ml, • • • , Ui-^-l ^Ui, • ' ' ) '^n) 

The energy of a multiparticle state is the sum of the one particle energies 

N 

E{ui, ■ ■ ■ ,un) = y^^ejuj). 

i=l 


In inhnite volume we normalize the in states as 


(5.2) 


(5.3) 


“(m'i, • • • • • • ,MAr)“ = {2 tt)^ 6 nm - u[) ■ --diuN - u'pf), (5.4) 

and the norm of a general state can be determined from (5.4) by (5.2). 

Let us consider the matrix elements of a local operator 0{t,x) between asymptotic 
states. The space-time dependence can be easily factored out 

°"*( m ' i ,--- ,u^M\0{x,t)\ui,--- = - , m ' m | 0 ( 0 , 0 ) | mi ,-•• 

where 

M N 

AF = ^e(u' ) - ^e(Mfc), (5.5) 

j=l k=l 

M N 

j=l k=l 

and we dehne the form factor of operator O as 

F'm,N ■ ,Un) = ■ ,u'm\ 0(0,0) |mi,- • • ,MAr)“ (5.6) 

^Although the states with non-ordered rapidities are not physical. 
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The form factors are a priori defined for ordered set of incoming and outgoing rapidities 
but can be analytically continued by (5.2). A form factor is a meromorphic function in all 
variables and each pole has a physical origin [34]. 

Suppose that the theory possesses crossing symmetry, i.e. a transformation which 
maps an outgoing particle with rapidity u to an incoming anti-particle with rapidity u. 
The crossing symmetry implies the crossing equation for the form factors which, in case of 
a single self-conjugated particle, reads as 

^M,N (^1) ■ ■ ■ ) ■ ■ ■ ) ('^1) ■ ■ ■ ) Ui, - ■ ■ , Un) (5-7) 

N k-1 

+ X] WM^k) n ^ {ui,Uk) («!,• • • ,Um-i\ui, ■■■ ,Uk,--- ,Un) 

k=l 1=1 

where the terms on the second line of (5.7) describe disconnected processes that occur if 
one of the incoming and outgoing particle has the same rapidity. The hat Uk denotes that 
Uk is missing from the list of rapidities. By using the crossing relation all form factors can 
be expressed in terms of elementary form factors 

F^(mi,--- ,un) = (0|C)(0,0)|rii, • • • ,un). (5.8) 

These elementary form factors satisfy several functional relations, called the form factor 
axioms. As these axioms are not relevant for the aims of this paper we will not list them 
here but rather refer to [34] for relativistic models and [19] for the AdS/CFT case. 

5.1.1 Diagonal form factors 

The diagonal form factor of the local operator O, dehned as 

‘“(tti,--- ,'UAr|C>(0,0)|rti, • • • (5.9) 

is singular due to the disconnected terms in the crossing relation (5.7). To avoid the 
singularities we regularize it by slightly shifting the outgoing rapidities. After crossing we 
get 


-Fei,-- - ,UN + eN,UN, - ■ ■ ,ui) = {0\O\ui + ,un + ei,ui, - ■ ■ ,nAr)“ (5.10) 

The diagonal limit, a —)• 0, of (5.10) is not well-defined. It was first noticed in [32] that the 
singular parts vanish in the limit when all —?■ 0, but the result depends on the direction 
of the limit. Its general structure can be written as 

(ui + ei, • • • ,UN + eN,UN, ■■■ ,ui) (5.11) 

N ^ N N N 

= ~ ^ ^ ^ ^ ^ ^ ("*^ 1 ; ■ ■ ■ ) Q1Q2 ■ ■ ■ + ■ ■ ■ 

i=l ^ ii = l 22 = 1 ^iV = l 

where is a completely symmetric tensor of rank N. The ellipsis denote terms which 

vanish in the Cj —?■ 0 limit. 
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There are two generally used regularization scheme in the literature. The first is the 
so-called symmetric evaluation when we set all Ci to be the same 

^2^® (ui, • • • ,un) = {ui + e,--- ,un + e,UN, ■■■ ,ui) ■ (5.12) 

The second scheme is called connected, in which the diagonal form factors are defined as 
the finite part of (5.11), i.e. the e-independent term, 

{ui,-■ ■ ,un) = N\ai2...N ■ (5.13) 

Both the symmetric and the connected diagonal form factors are completely symmetric in 
the rapidity variables u\, - ■ ■ Of course these two quantities are not independent and 

each can be expressed with use of the other [23]. 

5.2 Diagonal form factors in finite volume 

in this section we will summarize the results about the volume dependence of the diagonal 
form factors in all polynomial orders in the inverse of the volume, neglecting the exponen¬ 
tially small wrapping corrections, following [23, 35]. 

In finite volume L, the rapidities are quantized and a generic multiparticle state can 
be labeled by the Bethe quantum numbers - ,In})l- In finite volume we cannot 

order the particles by spatial separation in the remote past or future, as we did in the 
infinite volume case (5.2). In finite volume the states are completely symmetric under 
the exchange of particles and can be characterized by the set of quantum numbers. We 
adapt our notation to the conventions used in [23, 35] and order the quantum numbers 
in a monotonly decreasing sequence, Ii > • • • > /v The quantized rapidities with the 
quantum numbers {/i, • • • ,In} are solutions of the corresponding Bethe Ansatz equations. 
Similarly to (4.26) we define 

4>j({ui, • • • ,un}) = p{uj)L - i'^logS{uk,Uj ), (5.14) 

k = l 


and the quantization condition reads as 


■ ■ ■ ,un}) = 2-Klj, j = 1, ... ,N. 
These finite volume states are orthogonal to each other 


(5.15) 


• • • ,, • • • , In])l oc 5nm • • • ^In,Jn (5.16) 

and their normalization is a question of convention. 

One can change from the quantum number representation of states to the rapidity 
representation which gives the direct connection between the finite and infinite volume 

^ Apart from the free boson case all known S-matrix obey the property S{u,u) = —1 which is an effective 
Pauli exclusion. In this cases we have Ii > ■ ■ ■ > In ■ 
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states [35]. This change of variables involves the Jacobian, which is the density of N- 
particle states, defined as 


QNiui,--- ,un)l = detj^^\ui,-■ ■ ,un)l 


(5.17) 


= 


d^kjui, • • • ,un) 
dui 


kA = 1, 


,N. 


where we explicitly indicated the volume dependence of these quantities. The determinant 
(5.17) is closely related to the Gaudin norm of Bethe states (4.30)^ Then the relation 
between the infinite and finite volume states reads as 

1 


|{7i, • • • , In})l = 


,UN)LWi<jS{Ui,Uj) 


\ui, • • • , Un) 


(5.18) 


where the rapidities {ui} are the solutions of the Bethe Ansatz equations (5.14) correspond¬ 
ing to the quantum numbers {/i, • • • ,In}- This identification holds up to exponential cor¬ 
rections. The product of S'-matrices in the denominator ensures that the finite volume state 
is indeed symmetric under the exchange of particles. 

Defining the system in finite volume regularizes all the divergences appearing in the 
diagonal limit of form factors (5.11), thus the normalized finite volume diagonal matrix 
element 

- ,/Ar}|C>(0,0)|{Ii, • • • ,In})l 

is finite, completely well defined and does not depend on the normalization of states. How¬ 
ever, it is a challenging task to relate them to the infinite volume ones in the general case 
[30, 31]. The problem become considerably simpler if we neglect the exponentially small 
wrapping corrections. 

Up to wrapping, the finite volume A^-particle diagonal form factor (5.19) can be ex¬ 
pressed as a sum over the bipartite partitions of the full set {1,2, •• • ,N}, involving the 
infinite volume form factors and some kind of densities of states. As the diagonal form 
factors in infinite volume depend on the regularization scheme, this series is also scheme 
dependent. In case of the connected evaluation the relation reads as [23, 36] 

L{{h,-■ ■ ,In}\O{0,0)\{Ii,-■ ■ ,In})l 1 


• • • , 7Ar}|{/i, • • • , In})l 


Pn{{1, ■ 


,A^}) ^ 

’ ’ oC{i,...,Ar} 


{{Uk}k&a)PN (a) 
(5.20) 

where a denotes the complement of a in the full set. The functions appearing on the right 
hand side are exactly the connected diagonal form factors 


/^(ui, • • • ,ui) = 


0,C( 


,ui) 


(5.21) 


The functions are defined as the diagonal minor determinants of the A^-particle Jacobian 
(5.17), 

PN{a) = dei ■ ■ ,un)l , a C {1, • • • , iV}. (5.22) 

k,lGa ’ 

®The Gaudin norm itself is not physical as it depends on the conventions. However, in any convention, 
it is proportional to the Jacobian (5.17). 
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They can also be referred to as partial Gaudin norms. As special cases we have 


PNii'i-,-■ ■ = Qn{ui,-■ ■ ,un)l ; /OAf(0) = 1- (5.23) 


We want to emphasize that the function pN{oi) depend on all the N rapidities. The set of 
rapidities {ui} in the right hand side of (5.20) is the solution of the Bethe Ansatz equations 
(5.14) corresponding to the quantum numbers {/i,-- - ,In}- Thus, the explicit volume 
dependence is encoded only into the factors pN, the connected form factors depend on 
the volume only implicitly via the Bethe Ansatz equations. 

As the connected and symmetric diagonal form factors are not independent, we can 
express the finite volume matrix element in the symmetric regularization scheme. In this 
case the series take the form [35] 


L{{hr-- ,Jiv}|0(O,O)|{/i,-- - 

L{{h, • • • , lN}\{h, • • • , In})l 


1 

Piv({l,--- ,iV}) 


X] -^2^1 ({^fcifcea) P\a\ (a) • 


(5.24) 

Here again, the rapidities {ttj} are the solutions of the Bethe Ansatz equations (5.14) with 
the quantum numbers {Ii, • • • , In}- The p^a\ functions appearing in the sum are the |a|- 
particle densities of state (5.23,5.17) evaluated at the rapidities {ui}i^a- Note that, contrary 
to the connected expansion, they depend only on the rapidities labeled by the set a. The 
explicit volume dependence is carried only by the p functions. 


5.2.1 Form factor of densities of conserved charges 

An important special case of local operators is the density of a conserved quantity, 

Q= J{x, t)dx 

Jo 

where Q acts diagonally and additively on the multiparticle states. Its density therefore 
satisfies 

,/jv}|J(0,0)|{/i,-- - ,In})l _ 1 , 

l{{Ii,-■ ■ ,In}\{Ii,-■ ■ ,In})l ^ ^ ' 

where {ui} are the solutions of the Bethe Ansatz equations (5.14) corresponding to the 
quantum numbers {Ii, • • • , In}, and q(u) is the one-particle eigenvalue of the operator Q. 

A compact expression for the connected diagonal form factors of these densities was 
presented in [30, 36], however the proof was found recently [37]. The connected form factors 
can be cast into the form 

,Un)= ^ p\u„{l))p{Ua{l),U^{2))‘P{Uai2),Ua{3)) ■ ■ ■ P{Ua(N-l),Ua{N))q{Ua{N)) 

(5.25) 

where the summation runs over all the permutation of the set {1, • • • ,N}. Here, p' denotes 
the derivative of the momentum w.r.t the rapidity, 

/ ^ 

P{u) = -^Piu)- 
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For massive relativistic models we have p{u) = msinhtt and e{u) = mcoshu, so that 
p'{u) = e{u) and (5.25) reduces to the expression presented in [30, 36]. However, in the 
case of the Heisenberg XXX 1/2 spin chain an extra sign appears, p'{u) = —e{u). 

5.3 Conjecture for the symmetric structure constants 

Based on explicit calculations the authors of [1] conjectured that the Heavy-Heavy-Light 
symmetric structure constant in the su(2) sector of AA = 4 SYM is equal to a finite volume 
form factor and its volume dependence (up to wrapping corrections) at any coupling has 
the form (1.1). Let us suppose that the two, conjugated heavy operators correspond to 
a multiparticle state, labeled by the rapidities {ui}, in finite, but large volume L, such 
that the exponential corrections are negligible. Then the symmetric structure constant are 
conjectured to be the finite volume diagonal form factor of the vertex operator O of the 
light operator, 

Ca= ,un\0\ui,--- ,un)l (5.26) 

that can be expressed in terms of the infinite volume quantities (5.20) and (5.24). 

6 Matrix elements of spin operators 

In this section we study the diagonal matrix elements of spin operators of the Heisenberg 
spin chain using the ABA and the solution of QISP discussed in section 4. We show that, 
in general, the matrix elements can be written as linear combinations of a special kind 
of Slavnov determinant. In section 7, we show that this determinant has the structure 
conjectured in [1], namely it can be written as linear combinations of diagonal minors of 
Gaudin determinants (5.20). We call the procedure of expanding quantities in terms of 
diagonal minors of Gaudin determinants the finite volume expansion, as it captures all the 
finite volume dependence. We will discuss the case Ig = 1, 2 in detail and comment on the 
general ls> case. 

6.1 Form factors of length-2 operators 

We have shown in section 3 that all the length-2 diagonal matrix elements can be written 
as linear combinations of the following building blocks 

T""! = (u|oi(n)|u), = (u|o 2 (n)|u), i = l,...,3 (6.1) 

where the local operators oi(n) and 02 (n) are given in (3.9) and (3.10). According to (4.35), 
these matrix elements are proportional to the following quantities 


o^F^ = {u\A{en + i/2)\vL). 

(6.2) 

jrol ^ pAA ^ + i/2)A{en+l + i/2)|u). 

(6.3) 

j:ol ^ pBC _ + i/2)C{dn+l + i/2)\u), 


jol ^ pCB ^ + i/2)Bi9n+i + i/2)\n), 


j:o\ ^ pDD ^ + i/2)D{9n+i + i/2)|u). 
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In order to compute the building blocks (6.2) and (6.2), we act all the operators on the 
ket state |u). The action of A and D on a Bethe state is 


N 


A(?;)|u) = ^ Mn(?;) |{u,r;} \ {un}), 


(6.4) 


n=l 

N 


D{v)\u) = V ) + '^^niv) Ku,?;} \ 


n=l 


where Mfc(r;), Nfc(t>) are given by 


N 


Mt(o = jj 


Urj Uj % 


V -Un Un- Uj 

jjtn 


(6.5) 


N 


id{Un) Un-Uj+i 

^k{v) = - - z[[ 


Ur,. - V 




These relations can be derived from the Yangian algebra (4.7). From (6.4) we see that the 
action of the operators A and D on & Bethe state preserve the number of magnons. In 
addition to the original Bethe state |u), there is a sum of Bethe states |{u, v} \ {u^}) where 
one of the rapidities is replaced by the spectral parameter v of the operator. These are 
called the unwanted terms and are off-shell for generic v. On the other hand, they are not 
too far from the on-shell Bethe state |u) since most of the rapidities remain unchanged. 

The action of the C operator on the Bethe state is more involved 


N 


C'(^^)|u) = ^ |{u} \ {n„}) + ^ Kkn |{u, u} \ {Uk, Wn}) 


( 6 . 6 ) 


n=l 


k>n 


where 


ia{v)d{un) Uj — Un — i Uj — V + i 
Kn = - - --1- 

-1- 7/ .- — 7/„ 7/ .- — 71 


Ur,. - V 


(6.7) 


N 


ia{Un)d{v) -r-r Uj — Un + 1 Uj — V — i 


V — W 


n 


n . , 


Un — Ur. 


Ui — V 


Kfcn — 


In 

d‘(^Uf^'jCl(^Un') U}^ Un ^ p Uj U}^ % Uj Un ^ 


n 


{uk-v){un-v) Uk-Un Uj - Uk Uj - Un 

d{Un)a{uk^ Uk Un ^ T r Uj Uk i Uj Un ^ 
{Un-v){uk-v) Uk-Un Uj - Uk Uj - Un 

The coefficients and Kkn can be expressed in terms of and 

71 _ 7/ 7; _ 7/ 

K„(r;) =M„(7;)No(r;)-^ + Mo(r;)N„(t;)- 


+ 


V - Un +I 


Knk{v) = Mk{v)Kniv)—^ -+ Mn(u)NA;(t' 

Un-Uk + l 


V — Un — I 
Un Uj^ 


( 6 . 8 ) 


Un U}- % 
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if we define 


Mo = a(v) 


Qu('^ *) 

Qu(v) 


No = d(v) 


Qu(v + i) 
Qu(v) 


(6.9) 


From (6.6) it is clear that C reduces the number of magnons by one. For diagonal matrix 
elements, any C operator has to be accompanied by a i? operator in order to preserve 5^, 
and obtain non-vanishing results. For length-2 operators the only possibilities are B(u)C{v) 
and C{u)B{v). Both combinations preserve the number of magnons but will lead to a sum 
of unwanted terms with one or two magnons replaced by the spectral parameters of the 
operators. 


6.2 Form factors of length Ig operators 

In general, the action of Ig ABA operators on a Bethe state generates unwanted terms with 
at most Ig rapidities replaced by the spectral parameters of the operators. In particular, it 
is clear now that all the building blocks (6.2) and (6.3) can be written as the scalar products 
of the following three types 

(u|u), (u|{u,6»+} \ {wfc}), (u|{u,6»+,6»++J \ {rxj,ufc}). (6.10) 

where we have used the notation 0^ = -|- i/2. The first scalar product is the Gaudin 

norm (4.30). The finite volume dependence of other two determinants is the subject of 
discussion in section 7. 

The discussion of the previous subsection can be generalized to the diagonal matrix 
elements of the operators with Ig > 2. As before, any such matrix element can be spanned 
by some building blocks such as pBCA—A^ pD -D q£ number of the 

building blocks grows with the length of the operator. 

We act all the ABA operators on the ket state which give rise to the unwanted terms 

|{u,6'+,...,6'++^}\{ufci,...,nfc^}), M<lg. (6.11) 

Thus the diagonal matrix element of any length-/^ operator can be written as a linear 
combination of the following scalar products 

(u|{u,6»+,...,6'++^}\{ufc,,...,Ufc^J), M<lg. (6.12) 

The number of terms and the complexity of the coefficients will grow quickly with the 
increase of number of magnons and length of the operators, nevertheless the structure is 
robust. 


7 Finite volume expansion 

In this section, we analyze the structure of the scalar products (6.12) and show that any of 
them can be expanded in terms of diagonal minors of Gaudin norms. We call this procedure 
the finite volume expansion. 
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Above mentioned scalar products can be computed by the Slavnov determinant formula. 
In the Slavnov determinant (4.20), the denominator is a simple Cauchy determinant and 
can be computed readily. We therefore focus on the non-trivial numerator detj^ u^). 
Let us first consider the scalar product for the case of length-2 operators, (u|{u. On, 6n+i} \ 
{uj^Uk})- The determinant takes the following form 


i4>u ■ 

• • 

■ ^Ik ■ 

■ i 4>1N 

i(f)2i ■ 

• ^2j • 

■ ^2k ■ 

■ i 4>2N 

i4>Ni ■ 

• ^Nj ■ 

■ ^Nk ■ 

■ i4>NN 


(7.1) 


where we have defined i4>jk = ^{uj,Uk) and = Qik{ui,9n + il2). The procedure is 
straightforward: perform Laplace expansion with respect to the column or row that does not 
have any element of the form cfnn repeatedly, until one can not do it further. Note that after 
one Laplace expansion, we will obtain sub-determinants. We shall also perform the same 
procedure for all the sub-determinants until it terminates. This procedure will terminate 
when all the determinants in the expression take the form of diagonal minors (5.22) of 
Gaudin norm (5.17) 


PNiih, ■ ■ ■ ,im}) = (- 1 )”^ 



Therefore, the following expansion holds 

(u|{u, 6n, On+i} \ {uj,Uk}) = ^ J^{a) pn{ol), 

ckCA 


(7.2) 


(7.3) 


where A = {1,..., j,..., k, ..., N}^ and the summation runs over all possible subsets a of A. 
Here a is the complement of a in A. For an explicit and simple example, see Appendix A. 

Finally we need to justify why we call this procedure “finite volume expansion”. From 
a simple analysis below, it is clear that all the explicit L dependence are contained in the 
diagonal minors of the Gaudin norm. In the ABA approach, the diagonal matrix elements 
are given in terms of the following functions: the eigenvalue of the diagonal elements of the 
transfer matrix, a{u) and d{u), the products of functions f{u — v) and g{u — v) (6. 4,6. 6), 
and the matrix elements in the Slavnov determinant formula Lljk and <fjk. Under proper 
normalization, the functions a{u) and d{u) always appear in the expression as the ratio 
a{u)/d{u) = In fact, this kind of phase factor is either canceled by the same factors 

from the norm, or be replaced by products of scattering matrices using the Bethe Ansatz 
equations and they do not appear in the final expression. The products of f{u — v) and 
g{u — v) functions do not depend on L. The matrix element Lljk defined in (4.21) also has 

®Here j and k mean these two indices are absent. 
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no dependence on L. Finally, (pjk with j ^ k reads 


again, do not depend on L. The only dependence on L is hidden in the diagonal element 
4>nn- Recall that we have 


E 

m=l 


{Un - OmY + 1/4 


N 

E- 

Z = 1 
l^n 


/ / 'Pnl 


(7.5) 


In the homogeneous limit, where = 0 {m = 1, • • • , L), the first term of (7.5) becomes 
L/{u^ + 1/4) which depends linearly on L. When we perform the Laplace expansion, we 
carefully avoid expansion with this kind of terms and they only appear in the diagonal 
minor pjy^a). Therefore, the finite volume corrections are all contained in pi^{a). This is 
one part of the conjecture in [1] ' . 

We have shown in section 6 that any diagonal matrix element can be written as a linear 
combination of specific determinant with coefficients that do not depend explicitly on L. As 
we showed above, these determinants allow finite volume expansion, thus we can perform 
the finite volume expansion of any diagonal matrix element in the Heisenberg spin chain. As 
was shown in section 3, the diagonal matrix elements correspond to three-point functions 
of HHL type. Therefore we have shown that the structure of finite volume dependence of 
three-point functions conjectured in [1] is also valid at weak coupling at the leading order in 
the 5u(2) sector. In section 9, we will show that the structure also holds at one-loop level. 

This is only half of the story. In the conjecture [1], each coefficient J-{a) of p{a) is 
identified with the form factor of the same operator in infinite volume. In order to check 
this statement, it is desirable to have a formulation of the diagonal matrix elements of 
the Heisenberg spin chain directly in infinite volume. However, we are not aware of such 
a formulation, although it seems possible to do it in the framework of coordinate Bethe 
Ansatz. In principle, the infinite volume form factors for our case can also be obtained 
by first solving the Klose-McLoughlin axioms [19] and then take the weak coupling limit. 
However, no solution has been found up to now. Because of these reasons, we are not 
able to confirm that the coefficients we obtain from finite volume expansion are indeed the 
infinite volume form factors. 

It is still of interest to know the explicit form of coefficients from our finite volume 
expansion. These will be our predictions for the diagonal form factors in the infinite volume 
theory. We perform the finite volume expansion for all the diagonal matrix elements of 
length-1 and length-2 operators and extract the coefficients. The results exhibit a nice 
structure and will be presented in section 8. 

^There the authors used an equivalent description of the diagonal matrix element (5.26) in terms of the 
symmetric expansion (5.24), instead the connected one (5.20) that we used here. 
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8 Infinite volume form factors 


First, let us comment on the identity operator. As any multi-magnon diagonal matrix 
element of the identity operator equals to 1, matching it with the series (5.20), one can 
easily derive that all inhnite volume connected form factors vanish except from the vacuum 
expectation value, 

/(0) = 1 ; f{ui,...,UN) = 0 , A^>1. (8.1) 

We should also discuss separately the simple case of the vacuum expectation values of 
spin chain operators. In the series (5.20), the zero magnon diagonal matrix element only 
contains the vacuum expectation value of the given operator in the inhnite volume theory. 
So that, for the length-1 and length-2 operators one can easily hnd 

/°1 (0) = fol (0) = 1 , fol (0) = fol ^0) ^ ^ (8 2) 

It holds also in the general case. Let us take the operators E“^, (a, 6 = 1,2) as a basis 
on the local quantum space, and linearly extend it to Ig neighboring site. Then only one, 
among this 4^® basis element, has non-vanishing vacuum expectation value, namely the one 
containing at each site. 

In the rest of this section, we will perform the hnite volume expansion for the diagonal 
matrix elements of length-1 and length-2 operators. We will discuss a simple example, 
namely the case of length-1 operator with 2 magnons in detail and present the results for 
more complicated form factors. 


8.1 An example: length-1 operator with -2 magnons 

We consider the hnite volume diagonal matrix element for the operator oi(n) = E)j^ with 
two magnons 


T°^{UI,U2) 


{ui,U2\oi{n)\ui,U2) 

(ui,U2|ui,U2) 


It has the following structure in hnite volume 


(8.3) 


{ui,U2) = ^^(^1 (/32({1, 2}) (n2) P2({1}) + T' (^^l) P2{{‘i]) + («!, ^ 2 )) 

(8.4) 


where f°^ (u) is to be identihed with the connected diagonal form factor of oi in the inhnite 
volume theory (5.13). 

We proceed as described in the previous sections. Using the solution of QISP, we have 


^ tu{0n) {ui,U2\ui,U2) 


(8.5) 


where the denominator is the Gaudin norm (4.30), 


{UI,U2\UI,U2) 


( ^ 

Wa{uj)d{uj) 

Vi=i 


1 -F (ui - U 2 Y 

(ui - U 2 Y 


/^ 2 ({ 1 , 2 }). 


( 8 . 6 ) 
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From (6.4) 


{ui , U2\A{e^ ) |rii, ^ 2 ) = Mo (0+) (mi , M2 ki, U2) 

+ Mi((9+)(t6l,t62|u2,6'+) + M2((9+)(t6l,t62|ui,6'+). 


(8.7) 


We introduce some notations in order to simplify the expressions. Let us define 


Y\f=ia‘{vj)d{uj) 

uj-Vk+i 


V = 


( 8 . 8 ) 


so that 


(v|u) = det 0(Mj,Ufc). (8.9) 

jk 

By perform the hnite volume expansion for the three scalar products, we obtain 

{ui,U2\A{6^)\ui,U2) = -'^{ui,n2},{Ml,«2}^o(6*)))) P2({1, 2}) (8.10) 

- i'^{u„u,},{u„9i}^2{0tMu2,e+) P2({1}) 

- P2({2}) 

- </>i2 {^{u„u,Uu„eU^i^2, et)M,{e+) + 0+)M2(e)) 


Plugging (8.10) into (8.5) and comparing to the expansion (8.4), we obtain the expression 
for the various form factors in inhnite volume 


r{n,) 

r{u2) 


f°^{ui,U2) 


^ ^{ui,U2},{u2,ei} ^l(^n) r»^„, n+\ 

A/r 

®{«l,n2},{ni,«2} Mo(tln) 

} ^2{Gn ) Q/., fl+'i 
^{ui,U2},{ui,U2} Mo(vnj 


_^12_ 

^{u„u2Uu„u2}Moi0^) V i^^’-A,{u2,eU 

+‘*^W,.2},W,e+}^(«l>^n)M2(e)) 


Q{u2,B+)Mi{e+) 


( 8 . 11 ) 


Substituting the explicit expressions in (8.11) and, at the end, taking the homogeneous limit 
9n —t 0, we obtain very compact results for the inhnite volume connected form factors. 


n{nk)= 2 i-i/V ^ = ^’2 

% -I- J-/4 

" (nf + 1/4 + ul + 1 / 4 ) l + {ui-U 2 y 

8.2 Length-1 operator with N magnons 

We can perform the same calculation as in the previous subsection and extract the form 
factors with more magnons. The process becomes cumbersome for higher number of parti¬ 
cles. However, from the hrst few magnon cases, we are able to observe a nice pattern of the 
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connected form factors. The iV-magnon connected diagonal form factor for oi(n) is given 
as 


...,un) = £1 012 023+ permutations (8.13) 


where Sk is the energy of the magnon with rapidity Uk and (pjk can be seen as some “prop¬ 
agator” defined as 


£k = £{uk) = 


U\ 


+ 1/4’ 


01 fc — 


1 + {Uj - 


j^k. 


(8.14) 


The expression (8.13) can be represented by the diagrams in Fig. (2). Each node is labeled 
by a number from 1 to N. The leftmost node is associated with the energy of its label. 
The lines between two neighboring nodes are associated with a propagator. Multiplying 
the factors we obtain the value of the diagram. Summing over all the permutations of the 
labeling gives the result for infinite volume form factor f°^. The result for an N magnon 
state is thus a sum over N\ terms. 


031 012 024 

£3 •—o—o—o 

3 12 4 


Figure 2. Diagrammatic representation of one term in (8.13) with 4 magnons. 


The structure of is exactly the structure of the connected form factors of conserved 
charge densities (5.25). This is not surprising, since oi(n) = ^(I + a^) is indeed a length-1 
conserved charge density of the Heisenberg spin chain. The nice feature is that once we 
know the one particle eigenvalue q{u) of the charge, we can immediately write down the 
expression for the corresponding infinite volume form factors. We remark here that our 
result (8.13) is consistent with the determinant formula of [25]. 

8.3 Length-2 form factors of N magnons 

The calculation of infinite volume matrix elements can be performed following the same 
line as in section 8.1 but the process is more involved. Nevertheless, we again found some 
patterns for the various matrix elements which we present below. The structure for the 
length-2 operators can be encoded into diagrams similar to the one in Fig. (2). However, in 
this case we have two types of them, as are shown in Fig. (3). 

Each diagonal matrix element is given by two kinds of diagrams. This structure was 
tested for all length-2 operators up to 5 magnons and we conjecture that it holds for any 
number of magnons. The first kind is depicted in the blue region. We label the nodes 
by number from 1 to N. The leftmost node is associated with £{uj) while the rightmost 
node is associated with a function denoted by f^{uk), and it depends on the operator. The 
other kind, depicted in the green region, is more interesting. The leftmost and rightmost 
nodes are associated with Sj and e/ where £'{u) = In addition, for a given label of 

the nodes, one needs to sum over the diagrams which replaces one of the propagators by 
a “directed propagator”, = ip^{ui,Uj). The directed propagator is antisymmetric with 
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Figure 3. Different kinds of diagrams encoding length-2 form factors. 


respect to its arguments 00 = —0^ and its explicit form depends on the operator under 
consideration. 

In summarizing, the infinite volume diagonal matrix element of a length-2 operator, 
,U]sf), is characterized by two functions f^(u) and The result for N- 

magnon is given by 


f^{ui, ...,un) = (ei 0 i 2 --- 0 iv-i,AT + permutations) 


/N-l 


+ I ^ 61012 ... 0 j 0 j+i... 0 Ar-i,AT e)v TpermutationsJ 


. 2=1 


We list the data for o\ (3.9,3.10) is the following: 


r^u) = 2 
rhu) = - 
rhu) = - 
rHu) = 0 


u — f/2 
u + f/2 
u + f/2 
u — f/2 


v) = —{u — v){uv — l/4)0(u, v) 
0"2(u, v) = [u — v){u — i/2){v — i/2)4>{u, v) 

0"2(u, v) = {u — v){u -t- i/2){v + il2)(f){u, v) 
0"2(u,i;) = —{u — v){uv — l/4)0(u,u) 


(8.15) 


(8.16) 


Let us comment on this results (8.16). These data for the operators can be read off simply 
from the computation of 2 magnon case. Therefore, one should simply compute the 2- 
magnon matrix elements and perform the finite volume expansion to extract the data. 
Once the functions in (8.16) are known, we can write down any diagonal form factor of 
length-2 operators in the infinite volume. Any length-2 operator is a linear combination of 
the identity operator, oi and 02 , 

3 

On = bl + Co 0 i{n) + Co 0 i{n + 1) + Y,<^oUh) ( 8 . 17 ) 

2=1 

where b, cq and Cj (z = 1, 2, 3) are some numbers. Then the data of O is simply given by 

3 3 

f(u) = co + co + '^Cif^{u), = '^Ci'il;°2l^u,v), (8.18) 

2=1 2=1 
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and its vacuum expectation value is 

/^(0) = 6 + Co + Co + Cl- (8.19) 

For example, the operator o| is not independent 

02 (n) = I — oi(n) — oi(n + 1) + o^in). 

Note that by translational invariance, oi(n) gives the same result as oi(n + 1) when com¬ 
puting the form factors. It is easy to check that this resolution is consistent with (8.2,8.16). 
The diagonal matrix elements of the operators 02 and 02 are related by complex conjugation 
which is also manifest in (8.16). 

8.4 Examples of length-2 operators 

We compute two examples below. First, the matrix elements of the length-2 conserve charge 
density, which is the permutation operator or equivalently the Hamiltonian density 

Hfc,A:+i = Ifc,fe-i-i — Pfc.fc+i- We will see that the data for the permutation operator simplifies 
and the final result takes exactly the form predicted in (5.25). This is a non-trivial check 
of our functions (8.16). Another example is an HHL three-point function, with the light 
operator being the rotated BMN vacuum of length 4. 

Permutation operator The permutation operator is a length-2 operator of the 

Heisenberg spin chain. It can be written in terms of the operators O 2 with equal weights 

p/.,i=+i = 

i,j=l i=l 

According to (8.18), the data of permutation operator is given by 

f (u) = s{u) = ^|;^{u,v) = 0, /^(0) = 1. (8.21) 

The infinite volume form factor (8.15) with the entires (8.21) has the structure as a conserved 
charge should have (5.25), with the one particle eigenvalue of the corresponding charge being 
-e(u). 

A three-point function with Oa = TrZZZZ As another example, we compute a 
HHL three-point function with the light operator being Oa = O = Tr ZZZZ, which is the 
rotated BMN vacuum. The scalar field Z is defined as 

Z = Z + Z + i{X + X). (8.22) 

Following our strategy, we first map the field theory operator to the spin chain operator 
O using (3.4). Then we write the operator O in terms of linear combinations of the basis 
operators oi, O 2 (i = 1,..., 3). This enable us to write down the data for the operator O and 
thus the infinite volume form factor. We have® 

On = 4(1 — 2oi(n) — 2oi(n -|- 1) + 4o2(n) — 02 (n) — 02 (n)). (8.23) 

^Factor of 4 appeared due to the symmetry of the O 
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Using the representation (8.23) and also the data (8.16), we get for the operator O 

fiu) = =-6{u-v){4uv-l)(j){u,v). (8.24) 

The corresponding infinite volume form factors can be obtained from the general prescrip¬ 
tion (8.15) for the length-2 form factors. 


9 Matrix elements at one loop 


In this section, we generalize the above considerations from tree level to one loop. We 
show that the form factors at one loop can again be written in terms of a finite number 
of “building blocks”. These building blocks are matrix elements of the inhomogeneous 
Heisenberg XXX 1/2 spin chain with the inhomogeneities fixed to some specific values, called 
the BDS values, and can be written in terms of scalar products (6.12) for which one can 
perform the finite volume expansion, as at the tree level. 

There are several new features for three-point functions at higher loops. The dilatation 
operator in the su(2) sector is no longer the Hamiltonian of the Heisenberg XXX 1/2 spin 
chain, but becomes long-range interacting, called the BDS spin chain [38]. Therefore the 
two large operators correspond to the eigenvectors of the BDS spin chain. The BDS spin 
chain is related to a special inhomogeneous Heisenberg XXX;t /2 spin chain by a unitary 
transformation [11, 39]. Therefore its eigenstates can be obtained from the eigenstates of 
the inhomogeneous Heisenberg spin chain by performing the unitary transformation, namely 
|u)bds = S|u; 0®^®). The BDS values of the inhomogeneities are given by 

= 25sin^, k = l,---,L (9.1) 

where g is the coupling constant. The unitary operator S has been worked out in [11] up 
to order 


S = exp i 


k=i 




(9.2) 


where H^ = Ifc,A:-i-i ~ Pfc,fc-i-i and [Hj^ = [Hfc,Hfc_|_i]. The parameters and pk are related 
to the inhomogeneities by 

k k 

Uk = -^0j, pk = 2g‘^k-eki^k-'^0‘j, k = l,---,L. (9.3) 

i=i i=i 

It is obvious that Vk ^ g and Pk ^ g'^- 

Another new feature is that the quantum corrections manifest themselves as operator 
insertions at the splitting points [6, 7], as is shown in Fig. 4. The operator insertions take 
the form of the Hamiltonian density of the spin chain. In our set up we need insertions for 
the so (6) sector, since the light operator belongs to so (4), which is not closed at one loop. 

At one loop level, we need to take into account the aforementioned new features. 
In what follows, we first discuss the effect of the operator insertions and show that the 
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Figure 4. The quantum corrections are taken into account by operator insertions at the splitting 
points. At one-loop in the so(6) sector, the insertion takes the form of one-loop Hamiltonian density. 

three-point function under consideration can be reduced into the calculation of correlation 
functions of the BDS spin chain. Then we consider the effect of the S operator on the 
spin operators and reduce the correlation functions of BDS spin chain into the correlation 
functions of the inhomogeneous Heisenberg XXX 1/2 spin chain. As we discussed in the tree 
level case, the matrix elements of the inhomogeneous Heisenberg spin chain can be written 
in terms of scalar products (6.12), and we can perform the finite volume expansion. 

9.1 The effects of one-loop operator insertions 

For each spin chain state, there are two operator insertions at the two splitting points. We 
first discuss the effects of insertions for the “light” operator Oa- The one-loop insertion 
takes the form of the Hamiltonian density 

+ 2I;,;+1 - 2P;,;+1 (9.4) 

where and are the identity, permutation and trace operators. They act 

on the so(6) spin chain states as 

(9.5) 

6 

k=l 

At one loop level, the light operator should be an eigenvalue of the two-loop dilatation 
operator. The so(6) sector is closed only at one-loop so in principle one needs fields outside 
the so (6) sector, like fermionic fields, to construct the eigenstates of the two loop dilatation 
operator. However, when computing the three-point functions, the Feymann diagrams, 
involving fields apart from {X, Z, Y, X, Z,Y}, will not contribute at one loop order. Thus, 
the only new fields which we have to take into account at one loop are T and Y. It’s easy to 
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see from Fig. 5 that the light operator, in the presence of one-loop insertions, can again be 
written in terms of local spin chain operators , cr ^, due to the fact that the heavy operators 
are still can be expressed as su(2) spin chain states. Therefore, our considerations before 
can be generalized here. To summarize, the one-loop structure constants can be recast to 


zzzxzxxzzzzxz 



Figure 5. An example of mapping the light operator to the spin chain operator at one-loop in the 
presence of one-loop insertion. 

the calculation of matrix elements of the BDS spin chain BDs(u|0((T=^,cr^;5(2)|u) 

BDS- 

Now we consider the operator insertions for the heavy states. The effect of these 
insertions is increasing the length of the spin chain operator, as is shown in (6). 

This can be seen easily by noticing that 

2 

*J=1 


For example, we have the following 






Hull 


Hlul2 


(9.7) 


where the star stands for some strings of operators. 
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Figure 6. The effect of operator insertions for the heavy operators. They increase the length of 
the spin operator by 1. The red cross denotes the splitting point. 


9.2 The effects of the unitary S operator 

In this subsection, we discuss the action of unitary operator S on the spin operators. We 
are interested in the following quantity 


The S operator takes an exponential form S = expF, thus we have 


(9.8) 


S-ia+i(a±,(T^)S = di+i{a^,a^) - [F, d;+i(a±, + ^[F, [¥,di+i{a^,a^)]] +0[g^) 

(9.9) 

where we have truncated up to 0{g^) order. The action of S operator on the spin chain 
operator can be divided into two types. The hrst type is length preserving, it originates 
from the operators and [Hj^ that act within the range of the spin chain operator 0;+i, 
which gives rise to an operator with the same length, this is depicted in Fig. 7. 

L 


[HI 


; fTn T T r 

( 0;+l((T^,0-^) ^ 


. . i i ii iiii . . . . 

m • • • • • • m + l 


; TT T T T tr 


\ 


. . U ii i 


m +1 




Figure 7. The length preserving action generated by [H]*, on the spin chain vertex. 
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The other type of the action increases the length of the operator by 1 or 2, which are 
generated from the operators at the boundary of the spin chain operator. There are two 
kinds of length changing processes at one loop. One process is generated by a single or 
[H]fc, which is given in Fig. 8. 


VLiiTlT 


. . U i i i j. 



4 —♦— 


m + I 


;T T T T T T T t r, I 




.. m I T i- 


i • —• • 


m + / + 2 


Figure 8. The length changing action generated by [H]fc on the spin chain vertex. In this example, 
it increases the length of the spin chain operator by 2. 


The other process is generated by two H^’s, one example of which is given in Fig. 9. 



m + I 


t T T T T T T T T r | 

( j) 

m — 1 ••• ••• m + / + l 


Figure 9. The length changing action generated by two H^’s on both ends of the spin chain vertex. 
The length of the spin operator also increases by 2 in this example. 

From our analysis we see that the action of the S operator on the spin chain operators, 
in general, increases the length of the spin chain operator. Up to 0{g‘^) order, the length 
of the operator increases at most by 2. 

S-^di{a^,anS = 0;(a±,cT^) + + + 0(/) (9.10) 

This implies that in order to compute the form factor of length I operator for BDS spin 
chain, we need to compute the form factors of length 1 + 2, 1 + 1 and I operators in the 
inhomogeneous Heisenberg XXX 1/2 spin chain. 
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Once we write the three-point function in terms of matrix elements of the inhomoge¬ 
neous XXXi /2 spin chain, we can perform the finite volume expansion and organize the 
results in the form conjectured in [1], At one loop, the matrix element of Gaudin norm is 
modified. The equations (4.26), (4.27) and (4.29) are still valid, but the eigenvalues a{u) 
and d{u) are corrected 

L 

k=l 
L 

d{u) = n 

k=l 

where x{u) is the Zhukowsky map given by 

x(u) H-^ = u 

x[u) 

By replacing 

Uk + il2 _^ x{uk + i/2) 

Uk-il2 x{uk-il2) 

in (4.26) and expanding the result up to 0{g^) order, we obtain the Gaudin norm at one- 
loop, In fact, the replacement (9.13) gives the correct Gaudin norm up to wrapping 

orders[10, 11]. Our conclusion is that the structure conjectured in [1] is also valid at one 
loop level with respect to the one-loop Gaudin norm. 

Of course the coefficients or the infinite volume form factors at one-loop will be more 
complicated. It is an interesting question to see how the infinite volume form factors are 
deformed at one loop and whether it is possible to bootstrap to all loops. 

10 Conclusion 

In this paper we considered symmetric HHL correlators at weak coupling of the su(2) sector 
of AA = 4 SYM theory. Using the formalism of computing matrix element of XXX 1/2 spin 
chain (see e.g. [25]), we showed that at tree level, as well as at one loop, the finite volume 
dependence (up to the wrapping corrections) of this kind of correlator is given by the 
expression conjectured in [1]. Apart from giving the general arguments of the proof we 
computed the simplest non-trivial correlators at tree level, which correspond to the case 
when the length of the light operator is equal to 4. The structure of the coefficients of the 
finite volume expansion allowed us to conjecture their form for any number of excitations 
of the heavy operator. We showed that these coefficients can be expressed in terms of 
several functions of the rapidities defined by the light operator. Another aspect of the 
conjecture proposed in [1] suggests that the coefficients of the above mentioned finite volume 
expansion can be interpreted as appropriate infinite volume form factors. Therefore our 
explicit calculation provides us with concrete proposal for the infinite volume form factor. 



(9.13) 


(9.12) 


(9.11) 


- 32 





As a continuation of our work it would be interesting to see how the coefficients in 
the finite volume expansion are deformed at one loop. As an obstacle, one should mention 
that, as it follows from section 9, the complexity of the simplest one-loop computation 
is equivalent to a tree level computation with a light operator of the length equal to 6, 
and the calculation becomes much more involved. On the other hand, the recent result of 
[14] suggests an all-loop method for computing three-point correlators. And it would be 
interesting to merge their method with the approach developing in this paper and check 
the conjecture of [1] at all loop. In addition, it would give us all-loop prediction for the 
diagonal infinite volume form factors. 

Another interesting direction to pursue is to try to see whether the same structure of 
the finite volume dependence holds for symmetric HHL correlators of the other sectors of 
AA=4 SYM. 
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A An example of finite volume expansion 


In this appendix, we give an explicit example in order to illustrate how to perform the 
finite volume expansion of the special scalar products defined in section 7. The scalar 
product under consideration is ({rti, ri 2 , rts, rt 4 }|{tti, 0 , W 3 , M 4 }), where {wi, W 2 , ^* 3 ) ^ 4 } is a 
set of Bethe roots. Consider the numerator of the Slavnov determinant formula (4.20), 


({^tl,^t2,^^3,^^4}|{^^l,6',^t3,^i4}) OC 


(t>ll </>13 '/’14 

4)21 H 22 423 </>24 
431 ^32 433 </>34 
</>41 H 42 </>43 </>44 


(AT) 
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We first perform the Laplace expansion for the second column, which gives 



</'21 

4>23 

4’2A 


(pll 

(pl3 

4>1A 

f^l2 

4>31 

</>33 

</>34 

+ ^^22 

4>31 

4>33 

(p3A 


4>A1 

</>43 

(j)AA 


(pAl 

4>A3 

4>aa 



4>13 

4’ia 


(pll 

4>13 

4’ia 

f^32 

4‘21 

4>23 

4’2A 

+ f^42 

4>21 

4>23 

024 


4>AI 

4>A3 

(pAA 


4>31 

4>33 

034 


The Laplace expansion gives rise to 4 terms, which we shall denote Tj, z = 1, • • • ,4. For 
Ti, we do Laplace expansion by the first column 

Tl = — ni 2 {021/54({3, 4}) — </'3i((?^23/54({4:}) — </'24</’43) + </>41 (</'23</’34 “ </>24P4({3}))} 

= — ^124>21 /54({3, 4}) + ni2</>41</>23 P4({3}) + ril2</'31</'23 P4({4}) 

— ni2(</>31<?^24</>43 + (I>AI<^2a4>m) 

The second term already takes the form of diagonal minor of the Gaudin norm 

T 2 = 1^22 P4({1,3,4}) (A.3) 

For the third term, we perform Laplace expansion with respect to the second column 
T 3 = — ^32(t>23 P4({1, 4}) + n32<?^43</’24 P4({1}) + f^32</’13021 Pa{{^}) 

— 032((/>13(/)24</’41 + </’43</’14(?^21) 

For the last term, we perform Laplace expansion with respect to the last column 
T 4 = — VLi2(l>2A P4({1, 3}) + n42</’34</>23 Pa{{^}) + ^42</>14</>21 P 4 ({ 3 }) 

- 042((/>14(/)23</'31 + </'34</>13<?^21) 

Collecting terms from the above calculation, we obtain the finite volume expansion of the 
scalar product 

({tti, tt 2 , U 3 , Ui}\{ui,e, U 3 , U 4 }) (X (A.4) 

f ^22 /54({1, 3, 4}) — ^421^24 Pa{{^, 3}) — fl32</>23 Pa{{^, 4}) — Q.l24>21 /54({3, 4}) 

+ (^32</>43'^24 + ^A24‘3a4’23) P4({1}) + {^A24>1a4’21 + ^4l2<?^41</>24) /54({3}) 

+ {^I2(i>3l4’23 + ^432I?!>13</’2i) P4({4}) — f4l2((/'31</>24043 + 4>Ai4’234‘3a) 

— 032((/>13(/)24</'41 + </'43</>14021) “ f^42(</>14023</’31 + </>34</'13</’2l) 

In fact, it is not hard to convince ourselves that the similar expansion can be performed for 
general scalar products defined in section 7. For length-2 operators, we have the following 
finite volume expansion 

(u|{u,6ln,6ln+i} \ {uj,Uk}) = ^ Fap{a), A = {1,2, ...,N} (A.5) 

aCA 

In general, the terms of the expansion grows quickly with the number of excitations and 
the expansion coefficients might get quite involved. 
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